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The hypothesis is investigated, that the thermal structure of QCD phases at and near 
zero chemical potentials is determined by long range coherence, inducing the gauge boson 
pair condensate, and its thermal extension, representing a fundamental order parameter. 
A consistent model for thermal behaviour including interactions is derived in which the 
condensate does not produce any latent heat as it vanishes at the critical temperature 
inducing a second order phase transition with respect to energy density neglecting even- 
tual numerically small critical exponents. Localization and derealization of color fields 
are thus separated by a unique critical temperature. 
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1. Introduction 

The existence and nature of the QCD phase transition is theoretically accessible 
to a universal description underlying a thermal system with all its general and 
specific restrictions. At the moment key theoretical questions concerning the phase 
structure of thermal systems in QCD remain. Lattice QCD calculations lead to 
very clear predictions, establishing the lack of any phase transition at zero chemical 
potentials Th e va lidity of this method and its results is not uncontroversial, for 
example see Refs. 

There is a longstanding experimental effort to measure characteristic signatures 
of the phase transition in collisions of heavy ions at several facilities at AGS at 
BNL, SPS at CERN, RHIC at BNL, LHC at CERN and in the future FAIR at 
GSI. The present state of these studies yields significant results, well interpretable 
as a strongly interacting partonic matter formed in the earlier stages of the heavy 
ion collisions These results suggest the existence of a phase transition in the 
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course of the collision. Assuming the correctness of this conclusion, the detailed 
phase structure and order(s) of assocated transition(s) could not be determined 
conclusively so far . 

The experimentally measured reactions cover a large range of collision energies and 
chemical potentials in thermal analyses. In particular, with increasing energies the 
baryochemical potential in hadron-hadron or heavy ion collisions decreases. For a 
co mpa rison reducing these systems to a common baryochemical potential we refer 



In this paper we propose to give an "outline in principle" of the thermal phases of 
QCD at vanishing chemical potentials and their relation to the gauge boson pair 
condensate of QCD. We proceed to show that the two phases, separated by the 
critical temperature T cr at \i — can indeed be understood as a dynamical conse- 
quence of uniquely the bosonic pair condensate at nonzero quark masses, in contrast 
to supersymmetric models with no spontaneous breaking of supersymmetry. 
The two phases are represented approximate ely by a collection of noninteracting 
hadrons for < T < T cr and a collection of 8 gauge bosons and 4 flavours of quarks 
and antiquarks (u,d,s,c & anti) for T > T cr , in which interactions are modeled. 
We use here initially two collections of hadron resonances , denoted Ntype=65 
and Ntype=26 in the following. The Ntype=65 collection, on which most of the 
subsequent calculations of thermodynamic state functions rely, is an extension of 
the smaller Ntype=26 one. In the former ( Ntype=65 ) the 4 SU3fi meson nonets 
with J PCn = ++ , 1 ++ , 2 ++ ; l" 1 in the lowest mass p-wawe qq configuration 
are added, as well as the lowest mass open charm and anticharm s-wave pseudoscalar 
and vector mesons with composition qc and cq and baryons and antibaryons 
containing one charmed quark and their antiparticles together with overall s-wave 
composition qq'c and qq'c with q,q' = u,d,s are added, together with 3 
representative color gauge boson binaries (glueballs) with J PC = ++ , h , 2 ++ 
and two hybrid nonets of composition q g q' with J PCn = 1 h . 
The adopted approximation allows to illustrate the thermal phase properties in 
principle. 



2. Energy momentum density tensor and gauge boson condensate 

We begin with the local, symmetric and conserved energy momentum density ten- 
sor "& ^ v existing and renormalized as a basic consequence of QCD . The above 
properties imply exact Poincare invariance. 



to 




= i? 



} ( x ) ; d v 



= 



(1) 



The central defining quantities, of the existence and order of the phase transition, 
are pressure, energy density and entropy density for zero chemical potentials. The 
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local operator form of the trace anomaly in QCD reads: 

Sff (*)+ s o (x) 

5 o = (-2/3( 5 ) / 5 3 ) [ -!(:) F A (x) F>» A (a#)] s . d . (2) 
= -2b (s) F^ A (x)(:)] 

In Eq.[5]the subscript s.d. stands for " renormalization scale dependent", while bo = 
9/(167r 2 ). All quantities are defined to be renormalized and renormalization group 
invariant, except those with the subscript s.d. (3(g) denotes the (Callan-Symanzik-) 
rescaling function, where g is the coupling constant, j (:) F A v F M v A (:) is the 
local scalar density - to be called local gauge boson bilinear - composed bilinearly of 
field strength tensors F A v - the latter including multiplicatively the gauge coupling 
constant in their definition, relative to the perturbative field strength normalization. 
A denotes the color component within the adjoint representation of SU(3) C . (:)(:) 
stands for a suitable normal ordering, mj denotes the mass of quark flavour f : S ^ 
denotes the local scalar density of antiquark / and quark / flavors : 
Sjf(x) = (:) qfc S { x )l S f c { x ) ( : )i where s denotes the spin. 

The focus is to consider vacuum expected values of the local operators in Eq. [2j 
which by translation invariance are independent of the position x, suppressed in 
the following 

»T (ft| <V(*) M = - 2& o («| l(-)F A F^ A (x)(:) |ft) 

+ m f (ft | Sj f (x) |ft) [ ' 

( ft | (:)\F A Ft*" A (:) | ft ) = B 2 shall be called the gauge boson pair condensate , 
abbreviated by B 2 . 

( ft | S i f | ft ) induces spontaneous chiral symmetry breaking and is gen- 
erally called quark condensate. 

In connection with normal ordering ambiguities it is important to admit in the 
precise form of the energy momentum tensor a nontrivial vacuum expected value , 
which as a consequence of exact Poincare invariance must be of the form 

( ft | (x) | ft ) = -Tlnv Pvac 

{ Vfiv = dia 9 (!) -!> -1> _1 ) ; Pvac = ~Pvac } 

independent of x — > (4) 

(x) = <v (x) - (0]^ (x) |n> |ft) (n\ 

with 9 " A ■& p „ ( x ) = 0; ( ft | AiV (ac) | ft ) = 

In Eq. | ?si = I ft ) ( ft I denotes the projector on the ground state . 
From the two local, conserved tensors in Eq.Uonly /S.d fjLll (x) with vanishing vacuum 
expected value is acceptable as representing the conserved 4 momentum operators 
and their densities yielding the integral form 



% = td 3 x A^o (t,x) 



(5) 
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We use here throughout strictly thermal , 'extension in phase space' associated 
potentials, depending in subtle ways on vacuum condensates. To these potentials 
no vacuum associated spontaneous parameters like p vac = —p vac , defined in Eq. 
[4j contribute in a direct way, dominating in the limit T — > ; fi a = . 
From Eqs. [3] and [4] we obtain the relation and estimates 'pour fixer les idees' 



Pvac = B 2 + \ A 

32tt 2 



0.00302 GeV 4 
0.00658 GeV 4 



B 



0.125 GeV 4 121 (6) 
0.250 GeV 4 US 



A =-£/ m f (Q\ S ff |fi) 

~ / 1 ( \ m % + m\) = 0.00217 GeV 4 

3. Construction of a thermal model including interactions 

We follow the strategy layed out in Refs. ° taking into account the modifications 
described above, distinguishing two eventual phases 

1) the hadronic (h)-phase , with color localized within stable hadrons and selected 

hadron resonances corresponding to the two collections Ntype=65 C 26. 
Thermal potentials of the (h) -phase are approximate d b y those of free 
hadrons, neglecting decay widths, as described in Refs. EEl 

2) the quark-antiquark-gauge boson (qg)-phase, wherein thermal potentials are 

related but not equal to those of free quarks and antiquarks, restricted to 
the flavors u,d,s and c, and 8 gauge bosons pertaining to the gauge group 
SU3 C . Next we describe the modeling of interactions in the (qg)-phase, 
which deviates from noninteracting constituents assumed in Refs. l^JSI , 

We introduce for later use for the (qg)-phase , the Gibbs density g q°g and energy 
density g i°\g pertaining to noninteracting tricolored quarks , antiquarks with 
flavors 

u , d , s , c and eightfold colored gauge bosons 



/oo 
lEpdE 

S, = ( 2s P in <* q9 + 1) 



3 for q J 6 for q, q 
8 for g \ 16 for g 



(7) 



/3 = l/T;Z = Tlog[lTexp (-/?£)] 

A CO = - (d/d/3) gig (T) — T 2 (d/dT) g {0 q ] (T) 

In Eq.[7]the index a qg runs over the different constituents of the (qg) phase, while 
w a denotes the multiplicity beyond momentum phase space associated with the 
constituent a qg . The sign ( =p ) in the expression for I is - for bosons and + for 
fermions . 
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We choose the following masses for quark flavors u, d, s, c 

[ m u = 0.00525, m d = 0.00875 , g , 
m s = 0.175, m c = 1.27 ] GeV 
The absolute masses of t he u,d,s l ight flavors as well as their ratios m u : m d ■ 
m s = 3 : 5 : 100 , | H | 12 | 13 | ^ use( j nere pi av n0 decisive role in the present 
derivations, within generous ranges of ± 20 % . 

The inclusion of the charmed quark serves the purpose to check whether it has any 
significant influence on the thermal parameters in the region of T cr ~ 0.2 GeV , 
wich turns out to be in the few percent range . 

We proceed to modify the free quark antiquark gauge boson (qg-) parametrization 
of the Gibbs potential and the energy density, which for /i a = must obey the 
exact relation 

9q 9 = g qg (T = (3 - 1 ) ; - (d/dfi) g qg (T) = g e qg (T) (9) 

and g qg O g h , g e qg <H> g e h 
The Gibbs- and energy-densities in the hadron phase are constructed from the 
expressions analogous to the ones given in Eq. [7] , where the index a qg — > at, runs 
over a suitable choice of hadrons and hadron resonances as defined in Refs. with 
real masses and neglecting interactions among these states . 

The ensuing parametrization of interactions is understood as representing the phase 
structure in principle and not in numerical detail. The modeling of interactions in 
the qg-phase is performed setting two parameters k , A g , independent of tem- 
perature, as approximately parametrizing the interaction in the deconfined phase 
in a limited region of T > T cr 

g eqg (T;k) =k B f\ g {T) (10) 
gq g (T;k,Ag) = kgqg 0) (T)-Ag 

The parameter < k < 1 is taking into account the reduction of Gibbs 
density or pressure relative to the noninteracting (Stefan-Boltzmann) limit , noted 
in perturbative QCD calculations of thermal parameters for T ~ T cr of interest 
here j while the second parameter A g arises as integration constant from 
the differential equation ( Eq. |9] ) , which is clearly satisfied for arbitrary values of 
(k, Ag). 

We proceed in two steps to map out the structure of the phase transition 
I : the condition determining T cr k 

The equality of the energy den sities - in the hadron phase g e had for T < 
T cr as outlined in Refs. an d m the qg-phase as defined in Eq. [10] 
g e qg for T > T cr determine the critical temperature 

Qe had CO — Qe qg (T;k) o T = T cr (k) (11) 
The matching ( Eq. [IT] ) is further restricted to yield the value 

T cr - 0.2 GeV <->• k ~ 0.365 (12) 

in accordance with the estimate of one of us . 
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II : the condition avoiding singular behaviour of pressure gradient 

This condition implies 

9h (T cr ) = g qg {T cr ; k, Ag) O Ag = Ag {T cr (k)) (13) 

and determines Ag 

Ag ~ 0.0062 GeV 3 = 0.775 fm ~ 3 ( 14 ) 



4. Results and discussion 

We begin with a disussion with a comparison of the two colletions of hadron reso- 
nances Ntype=65 and Ntype=26 whence thermal state functions are extrapolates 
beyond the temperature range 130 MeV < T < 170 McV obtained in studies 
of chemical equilibrium of hadron yield ratios. 

It is expected that approaching a phase transition will require interactions on the 
hadronic side as well as on the quark-gluon side. The mere extension of the collection 
of hadron resonances also beyond Ntype=65 here , without introduction of interac- 
tions among the hadron resonances is unrealistic. If there is no phase transition at 
all, this does not change the mentioned expectation. 

This is illustrated in Figs. 1 and 2 , in which the two noninteracting collections of 
hadron resonances, Ntype=65 and 26, are compared with respect to the quantities 
dscale = ( g e — 3 p ) / T 4 and pressure / T 4 respectively . 




T [GeV] 



Fig. 1. The thermal quantity representing the trace anomaly - dscale (T) = ( g e — 3p) / T 4 
is shown for 0.1 GeV < T < 0.7 GeV . The comparison of HRG collections Ntype=65 and 
26 shows the sensitive temperature regions beyond temperatures where chemical freeze out takes 
place - T ~ 150 - 170 MeV for Pb - Pb collisions at SPS and Au - Au collisions at RHIC. 
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Fig. 2. Pressure / T 4 versus T in the range 0.1 GeV < T < 0.4 GeV . for the t wo collections 



Ntype=65 and 26 is compared with the results of lattice QCD by S. Borsanyi et al. 



28 



As a consequence of the fair agreement of the pressure / T 4 as obtained from the 
noninteracting hadron res onanc e collection Ntype=65 with the lattice simulation 
results of S. Borsanyi et al. BUS] ^ ag s j 10wn m pjg 2 , we choose this collection in all 
subsequent calculations representing the hadron phase for T < T cr ~ 200 MeV . 
We mention the fact, that collections of hadron resonances used in recent analyses of 
especially ratios od hadron yields as observed in Au-Au collisions at RHIC include 
all resonances reported by the PDGQ^up to a mass of 2 GeV or 2.5 GeV , which 
are larger than the one defined Ntype=65 here. In Ref. ^ the agreement up to 
T ~ 170 MeV of all thermal state functions with the HRG quantities of these 
large collections is very satisfactory . 

The restriction to noninteracting hadron resonances nevertheless has a limited va- 
lidity particularly at temperatures above chemical freeze out as T ~ 200 MeV . 
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The result of steps I and II above are illustrated in Figs. 3-7 below . For a 
comparison with thermodynamic analyses in lattice QCD we refer to Refs. 

and The existence of two crossing points, one each for energy density and 
pressure, coincident in temperature at acceptable values of all parameters near 
T cr ~ 0.2 GeV is nontrivial , shown in Fig. 3 below 




-0.01 ' ' 1 ' 1 1 1 

0.1 0.12 0.14 0.16 0.18 0.2 0.22 

T [ GeV ] 

Fig. 3. The two thermal quantities energy density q e and Gibbs density g = p / T in units of 
GeV 4 and GeV 4 respectively are decomposed into the hadronic (had) ( valid for T < T cr ) 
and quark-gluon (qg) parts ( valid for T > T cr ) each : q e had : blue curve ; g e qg : purple 
curve and g had : red curve g qg '■ green curve . The two small circles mark the two equalities 
8 e had = 8 e qg and g had = g qg coincident in temperature T = T c ( = T cr ) . 



In Figs. 4-6 the modifications of g e qg /T 4 , p qg/T 4 and dscale = 
( Q e qg — 3 P qg ) / T* introduce a subleading critical exponent v ; v — 0.975 
used here throughout, in the vicinity of T = T cr allowing the free quark-gluon limits 
to be reached for T — > oo . The modified quantities are defined as 



Qeqg \ = / fmod ( v , T / T c ) Q e qg + fmodl ( v , T / T c ) p 



I f fmod ( v , T I T c ) Q e q 
J \ fmod ( v , T I T c ) p qg 



fmod ( v , T I T cr ) = 1 + | 1 - T cr I T \ 2u ( 1 / k - 1 ) 
fmodl {v 1 T/T c )=T{d/dT) fmod ( v , T / T c ) 



(15) 



Assuming that all thermodynamic potentials are analytic functions of T in a suffi- 
ciently large neighbourhood of T cr , with the transition arising from appropriate 
conditional choices below and above T cr , as is the case for the parametrizations 
adopted here, the orders of the so deduced phase transition (for \i a — ) 
become 
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order with respect to quantity 



2. f-> energy density 

3. -f-> pressure 

1. f-> sound velocity (16) 

We note the behaviour of the square velocity of sound. Steps I and II imply that 
the sound velocity (square) exhibits a 1. order transition at T = T cr . 
Let the squares of the sound velocity in the hadron phase ( qg-phase ), always for 
vanishing chemical potentials, be denoted respectively 



Phad (qg) 

J s had (qg) 

Qe had (qg) 



<hadiaa^ CO = ! ' = d/dT (17) 



It follows that the sound velocity (square) jumps from a lower value in the hadron 
phase to a higher value in the qg-phase . 
For the qg-phase we have the relations 

P q9 = T (kg qg 0) -&g) , TgM=pM (T) 

■ (0) (0) (rp\ entropy density 

V qg Qs \J- ) '• for nonintcracting 

qg-phasc (1^) 



vi 



.,2 _ * " - - I .,2 \(°) 



q9 ~TJ» 99 ' kd {0) 

t Qe qg n fe qg 

All quantities with superscript refer to the noninteracting qg-phase . The 

introduction of interactions leads to a reduction of the qg-sound velocity square 
proportional to A g . This is however a small correction for all T > T cr as 
shown in Fig. 4, also illustrating the 1. order transition of « s 2 .In the qg-phase 
the sound velocity square approaches the limiting value ^ and is just after the 
transition already quite near this value. The behaviour of « s 2 is shown in Fig. 7 . 
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l Nlype-65 

: i / 


/ 


--^ i— 


qg-limrt 


: \ / 




> / 


rho {e qg} / 1*^4} mod 

nu =0.975 - 


: \ / 




: 






rho {e qg) ! T*-{4) no mod 


rtio Jehgf/T"{4} 


T [GeV] 



0.15 0.2 0.25 0.3 



Fig. 4. The ( dominantly ) second order nature is displayed of the transition - unmodi- 
fied and modified on the quark-gluon-side by subleading critical exponents — of the quantity 
Q e ( T ) / T 4 , piecewise described by g e h 9 (_ 6 5) / T 4 for T < T cr from the HRG 
side and p e qg no mod (mod) / T 4 for T > T ar from the quark-gluon side as defined in Eq. 
I15I The three curves are represented with dashed lines outside their range of validity. 



Ntype 



/ 



qg -limit 



/ 



/ 



P_) qg) / T»{4] mod 
nu =0.975 




p_(hg}/T-{4] 



T [GeV] 



Fig. 5. The ( dominantly ) third order nature is displayed of the transition under the same 
conditions as underlying Fig. 4 for the quantity p / T 4 , piecewise described by p e hg(— 65) / T 4 
for T < T cr from the HRG side and p qg no mod ( mo j\ / T 4 for T > T cr from the quark- 
gluon side as defined in Eq. I15I . The three curves are represented with dashed lines outside their 
range of validity. 
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5 | 1 1 1 1 1 1 5 

\ / 

4-5 - \ J - 4.5 




Fig. 6. The ( dominantly ) second order nature of the transition is displayed with respect to the 
therma l qu antity — dscale - forming the shape of an 'indian tent' . The same quantity as obtained 
in Ref. from lattice simulation of QCD under the same thermal conditions is also plotted for 
comparison . The quantity dscale v qg mod is obtained from Eg. 1151 . here with v = 0.975 , the 
same as in Figs. 4 - 5 . The three curves forming the 'indian tent' are represented with dashed 
lines outside their range of validity. 
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Fig. 7. The ( dominantly ) first order nature is displayed of the transition with respect to the 
square of the velocity of sound , under the same conditions as for Figs. 4 - 6 , piecewise described 
by the quantities v ? ; ( hq = h ) for T < T cr from the HRG side and v 2 . for 

« ^ ng > \ & / — qg no mod 



curves for v ? 
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are represented with dashed lines outside their range of validity. 
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The derivations concerning the phase structure of QCD at vanishing chemical po- 
tentials , laid out in the present discussion , follow within the hypothesis , that the 
gauge boson pair condensate is at the origin of stable embeddings of quark flavors 
light and heavy - into the ( long range ) dynamics controled primarily by gauge 
boson self-interactions . By the same hypothesis - this embedding depends in a 
minimal way on quark masses , including their chiral as well as anti-chiral limits : 
m X i>'" m XM ; m h 1 , ■ ■ ■ m h N ->• oo . 

Extensive lattice simulations of the thermodynamic state in particular at zero chem- 
ical potential have been performed in recent years. With realistic quark masses for 
the three light flavours u,d,s no phase transition but a cross over region is found. 
For a review we refer to Ref. [13]. 

For pure gauge theories with gauge group SU N c , N c > 2 a 1st order phase transition 
and for N c — 2 a second order transition is found [14,15]. 

The latter phase transition agrees with the results pertaining to the 3-dimensional 
Ising model. 



The argument for - at least - a second order phase transition with 
respect to energy density 

In this subsection a short outline is given of how spontaneous parameters at T = 
and T > are obtained in general renormalizable local field theories. 
We define first in physical space-time, the effective action T for the composite clas- 
sical field ipd O <p associated with the local gauge boson bilinear, introduced in Eq. 

m 

exp iW = (ft|exp [iJd A x{C + Jtp) (x)j |fi) c 

= l(0*£ (?) FftV A (z)(0 ; J(x) ■ ^urce 

in the presence of a classical local source J . c in Eq. [T2] denotes the connected part 
of the generating functionsal W. The effective action T((p c ) is obtained from the 
generating functional W through functional Legendre transform 

<p cl (x) = 6/6J(x) W(J) ->• 

T((f c i) = J d 4 x ((p c i(x)J(x)) - W ; ipd = (fd(J) (20) 
8/5cp cl (x) T = J(x) ; <Pd(x)\j^ - B 2 , (Eqs.EH) 
T(ip c i) = / d^x v((p c i) 



The (nonlocal) density functional v defined in Eq. [20] determines the spontaneous 
ground state expected value of ip c i in the limit of vanishing sources as absolute 
minimum. 

It is essential to safeguard appropriate boundary and/or integrability conditions for 
the sources J{x) which translate into corresponding conditions for the Legendre 
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transforms <p c i(x), e.g. 

lin^^oo ./(a;) = Jd 4 y \J\ 2 (y) < oo (21) 

as described in Ref. ^1 The functional W , T in Eqs. [121 HOI can be extended to 
Euclidean space. The absolute thermodynamic limit at zero temperature is reached 
from the boundary conditions in Eq. [5TJ for spontaneous vacuum p aram eters al- 
lowing sources and classical fields to approach also constant values whereas 
thermal environment at finite temperature T corresponds to a finite Euclidean time 
extension 

A t E = = l/T (22) 

In the finite temperature thermal limit new boundary conditions at the bounding 
times tg — , /3 have to be set, periodic for gauge invariant bosonic sources and 
classical fields, as considered above relative to the local operator ip(x) defined in Eq. 
IIUI Hence we face the two thermodynamic limits, treated here for just the associated 
quantities 



|(:)i^ (x) F»» A (!)(:) = <p(x) -> ( cp cl (x),J(x) ) 
J(x) = (l/gi)-(l/gi(x)) = -A^-(x) 



(23) 



g 

The external source J{x) represents a space time dependent coupling constant, 
before thermodynamic limits are taken. 

Absolute thermodynamic ( T = ) - and thermal finite temperature ( T finite ) 
limits yield the relations 

J(x) _^Je[-J»,J»]^0 

v(tp c [(x)) T v(ip c i,T) 



T 




(24) 



d Vcl v{(pd,T) = J -)• -> ipd=ip c i(T) 



Performing the sequence of limits and letting the (constant) values of the source J 
vary in a suitable interval [— J*, J*] before relaxing it to zero, as shown in Eq. 1241 
defines the hysteresis line of the spontaneous parameter (p c i(J — > 0,T). The latter 
emerges as thermal average for T finite and as vacuum expected value for T = 
(Eqs. HUT]). 

In the T finite case we have the choice to include chemical potentials , one each for 
conserved quark flavor neutral currents, which we set to zero here. 
The phase transition underlying the present discussion is associated - by hypothesis 
- with the vanishinging in a singular way at finite critical temperature, T — T cr > 
of the quantity tp c i (T) defined in Eq. [24] 



Pd(T) for T -^T c 



(25) 
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Classical configurations leading to this behav iour in herit a 'Watt-less' nature from 
beeing particular to the ground state of QCD (22]23] anc j thus should not generate a 
step-like behaviour of thermal energy density. This is borne out in the subtraction 
of the ground state projection of the associated energy momentum density tensor 
-d al/ as shown in Eqs. 14151 

This however does not imply that the first derivative of energy density changes by a 
finite amount through the transition, i.e. its corresponding genuinely second order 
nature. The singularity could well be characterized by critical exponents without 
relation to a definite step in a given derivative - of thermal energy density. 

The mechanism of gauge boson pair condens ation d iscu ssed here is a new analysis 
taking its roots in material presented in Refs. ^1 and in chronological order. 
It is centered on the thermal average of the field strength bilinear local density 
operator defined in Eq. [2] 

(X) F»» A (*)(:) (26) 

The nonzero vacuum expected value ( Eq. [3] ) signals a connection between the 
localization of color and the structure of Bogoliubov transformations as appropriate 
for gauge fields at long range. 

For QCD, actual calculations of effective potentials at large and intermediary range 
are not accessible to perturbation theory. Hence analytic control is lacking at 
present. 

5. Conclusions 

Our modeling of strong interactions in the qg-phase, in the vicinity but above the 
critical temperature T cr ~ 0.2 Mev, as described between Eqs. 1101 andll8| is sup- 
ported b y obs ervations of strong coupling in central Au-Au collisions at 200 GeV 
at RHIC I312ZI. The reduction of pressure and energy density relative to noninter- 
acting (anti-) quarks and gauge bosons then allows for the phase transition to be 
(essentially) of second order with respect to energy density, as demonstrated in the 
related thermal quantities shown in Figs. 2 and 3. The square of the velocity of 
sound reveals most directly the order of the transition through its discontinuity as 
shown in Fig. 4. 

The upcoming lower cm. energy scan at RHIC, new results expected from LHC 
concerning hadronic physics and more extended studies in lattice QCD promise to 
clarify its phase structure, eventually accompanied by theoretical insights. 
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